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We study cosmic shear statistics using the ray-tracing simulation combined with a set of 
large A'^-body simulations. In this contribution, we first describe our method. Then we show 
some selected results especially focusing on effects of the deflection of light rays and the lens- 
lens coupling which are neglected in making the theoretical predictions of the cosmic shear 
statistics. 
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1 Introduction 

The cosmic shear statistics have been known as a powerful tool for probing the large-scale 
structure formation as well as for placing a constraint on values of the cosmological model 
El. Recently, four independent groups have reported detections of cosmic shear variance oQffl. 
Although those detections were made with relatively small fields which limit their accuracy, on 
going wide field cosmic shear surveys will provide a precious measurement of not only the cosmic 
shear variance but also higher order statistics such like the skewness of the lensing convergence. 

Since the pioneering work by GunnQ, there has been a great progress in the theoretical study 
of the cosmic shear statistics Q. The analytical formulae for making the theoretical prediction 
of the cosmic shear statistics are based on the perturbation theory of the cosmic density field 
combined with the nonlinear clustering ansatz. The accuracy and limitations of the theoretical 
predictions should be tested against numerical simulations that is one purpose of numerical 
studies &□. 

In this contribution, we summarize, briefly, the methods and selected res ults o f our research 
project on the numerical study of the cosmic shear, details are presented inll2lO. The project 
aims (1) to test the analytical predictions against the numerical simulation, (2) to examine the 
higher order statistics of the cosmic shear, (3) to simulate a cosniic shear observation to examine 
possible systematic effects caused by, e.g., the source clustering EJ'EJ. 



Table 1: Parameters in Three Cluster normalized CDM Models 
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Figure 1: Tiling configuration of simulation boxes. Dashed lines shows the comoving angular diameter distance 

of ±2.5 degrees. 



2 Models and methods 



We consider three cluster normalized cold dark matter (CDM) models, parameters in the models 
are summarized in Table Q. A^-body simulations were performed with a vectorized particle-mesh 
code. They use 256^ x 512 particles and the same number of force mesh in a periodic rectangular 
comoving box and use light-cone output ll3o. In order to generate the density field from z = 
to z ~ 3, we performed 11, 12 and 13 independent simulations for SCDM, OCDM and ACDM 
model, respectively. We adopted the tiling configuration of the boxesH, i.e., the box size of each 
realization is chosen so that we have a field of view of 5 x 5 square degrees (see Figure |l]) . 

Light ray trajectories are followed throiigh the density field generated by A^-body simulations 
adopting the multiple-lens plane algorithm t3. The lens planes are located at intervals of 80h~^. 
The initial ray directions are set on 512^ grids with the grid spacing of 5°/512 ~ 0.59 arcmin. For 
each ray, positions of the ray on each lens plane are computed, and then the lensing magnification 
matrix, Mij, is computed at the ray position on each plane. The lensing convergence, shear and 
net rotation are expressed by k = (Mn + M22)/2, 71 = (Mn - M22)/2, 72 = {Mu + M2i)/2, 
and u = (M12 — M2i)/2, respectively. We performed 40 realizations for each model changing the 
underlying density field (i.e., making random shifts of boxes to x and y directions (perpendicular 
to the line-of-sight) using the periodic boundary condition in A^-body simulations). 

Figure |2| shows the lensing shear map overlaid on the convergence map for ACDM model, 
the sources are assumed to be at a single redshift of Zg = 1. The filed is 5-degree on a side 
and contains 512^ line-of-sights. The angular resolution is limited by the spatial resolution of 
A^-body simulation. We found that the effective resolution is about 2 arcmin for the source 
redshift of Zs = 1 and is slightly better (worse) for the higher (lower) redshift. 




Figure 2; The lensing shear map overlaid on the convergence map for ACDM model, the sources are assumed to 
be at a single redshift of Zs — 1. The filed is 5-degree on a side and contains 512^ line-of-sights. 
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Figure 4: Top panel: The skewness parameter S3 as a 
function of the filter scale of top-hat window. Bottom 
panel: The difference in S3 between measured from 
K field of full ray-tracing and that of approximated 
ray-tracing. 



3 Results and discussion 



Figure y shows the root-mean-square of the top-hat filtered lensing convergence measured from 
the ray-tracing simulations compared with the nonlinear predictionEJ. The dumping on smaller 
scales comes from the fact that the lack of the power on scales smaller than the resolution in 
A'^-body simulation, while that on larger scales reflects the finite field effect, i.e., a lack of powers 
on scales larger than A^-body simulation box. One may find very good agreement between the 
measurements and the nonlinear predictions. 

The skewness parameter 5*3 defined by 6*3 = {k^) / {k^)"^ is known as a powerful probe of 
density parameter, r^mEJ, though, it is difficult to make a precious prediction of its value because 
a description of the nonlinear evolution of the density bispectrum is required. The upper panel 
of Figure ^ shows 6*3 measured from the results of ray-tracing simulations compared with the 
prediction based on the quasi-linear theory of density perturbation. It is clearly shown in Figure 
4 that the nonlinearity of the evolution of density field is very important even at the scale ~ 10 
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In addition to the usual full ray-tracing simulations, we performed approximated ray-tracing 
simulations to examine the effects of the defiections of the light rays and the lens-lens coupling 
Ii3 which are neglected in making the theoretical predictions. The procedure of the approximated 
ray-tracing simulations is same as that of full ray-tracing except for artificial by-passing of both 
the defiection of light rays and all lens-lens couplings. 

It is found that the difference in the variance of the lensing convergence (and of shear) 
between measured from full and approxim,ated ray-tracing is very small, \/S.RM S / RM S\ < 0.01. 
Therefore the the defiections of the light rays and the lens-lens coupling can be neglected safely 
for making the theoretical predictions of the variance. For the convergence skewness, we found 
that the difference is not significant, |AS'3//S'3| < 0.05 for three cosmological models. These 

americal results are compared with the theoretical predictions and a good agreement is found 
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